fc-string tensions in SU(A^) gauge theories 
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In the context of four-dimensional SU(A'') gauge theories, we study the spectrum of the confining 
strings. We compute, for the SU(6) gauge theory formulated on a lattice, the three independent 
string tensions related to sources with Zn charge fc = 1, 2, 3, using Monte Carlo simulations. Our 
results, whose uncertainty is approximately 2% for k = 2 and 4% for fc = 3, are consistent with the 
sine formula a^/o- = sin fc.^/ sin for the ratio between at and the standard string tension a, and 
show deviations from the Casimir scaling. 

The sine formula is known to emerge in supersymmetric S\J{N) gauge theories and in M-theory. We 
comment on an analogous behavior exhibited by two-dimensional S\J{N) x S\J{N) chiral models. 

PACS Numbers: 11.15.-q, 12.38.Aw, 12.38. Gc, ll.15.Ha 



Quantum Chromodynaniics is a nonabelian gauge the- 
ory based on the gauge group SU(3). The mechanisms 
underlying many of its fundamental properties, such as 
confinement, chiral symmetry, topological effects and the 
axial anomaly, are under active investigation; they are 
being studied by different approaches, including numer- 
ical simulations of the theory formulated on the lattice, 
several models of the vacuum, as well as some recent 
proposals derived from M-theory and AdS/CFT. Many 
features of QCD can be better understood by extending 
the study to SU(A^) gauge theories with N larger than 
three and in particular by examining the large- iV limit. 

Four-dimensional gauge theories exhibit confinement, 
i.e. static sources in the fundamental representation de- 
velop a linear potential characterized by a string tension 
a. As pointed out in many studies, it is important to 
investigate the behavior of the system in the presence of 
static sources in representations higher than the funda- 
mental one. This may provide useful hints on the mecha- 
nism responsible for confinement, helping to identify the 
most appropriate models of the QCD vacuum and to se- 
lect among the various confinement hypotheses. Among 
the latter, the so-called Casimir scaling hypothesis for 
the potential between heavy-quark sources in different 
representations has attracted much interest (see e.g. the 
recent publications 

SU(A^) gauge theories confine by means of chromoelec- 
tric flux tubes carrying charge in the center of the 
gauge group. A chromoelectric source of charge k with 
respect to Z^ is confined by a fc-string with string ten- 
sion (7k (cti = ct is the string tension related to the fun- 
damental representation). If ak < k a, then a string with 
charge k is stable against decay to k strings of charge 
one. Charge conjugation implies Ufc = un-u- Therefore 
SU{Z) has only one independent string tension determin- 
ing the large distance behavior of the potential for fc 7^ 0. 
One must consider larger values of N to look for distinct 
fc-strings. 

As pointed out in Ref. it is interesting to compare 



the fc-string tension ratios 



R{k,N) = — 
a 



(1) 



in different theories. The idea is that such ratios may 
reveal a universal behavior within a large class of models 
characterized by SU(A^) symmetry, such as SU(A^) gauge 
theories and their supersymmetric extensions, ft has 
been noted that stable /c-strings are related to the totally 
antisymmetric representations of rank fc, and that in var- 
ious realizations of supersymmetric SU(iV) gauge theo- 
ries i?(fc, A^) satisfies the sine formula i?(fc, A^) — S{k, N) 
where 



S{k,N) 



sin{kTT/N) 
sin(7r/Af) 



(2) 



R{k, N) has been computed for the JV = 2 supersym- 
metric SU(A^) gauge theory softly broken to M ^ 1 
obtaining Eq. (^. The same result is found also in the 
context of M-theory, and extended to the case of large 
breaking of the M = 2 supersymmetric theory [^j. An 
interesting question is whether the sine formula holds 
in nonsupersymmetric SU(A'^) gauge theories. The M- 
theory approach to nonsupersymmetric QCD, although 
it is still at a rather speculative stage, suggests that it 
may be so [^,||. However, as discussed in Refs. cor- 
rections from various sources cannot be excluded, so that 
this prediction cannot be considered robust. 

Another interesting and suggestive hypothesis is that 
the fc-string tension ratio satisfies the so-called Casimir 
scaling law |l0|, i.e. R{k,N) = C{k,N) where 



C{k,N) 



k{N -k) 
N ~1 



(3) 



is the ratio between the values of the quadratic Casimir 
operators in the rank-fc antisymmetric and in the funda- 
mental representations. The Casimir ratio is satisfied on 
the one hand by the strong-coupling limit of the lattice 
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Hamiltonian formulation of SU(iV) gauge theories [0, 
and on the other hand by the smaU-distance behavior 
of the potential between two static charges in different 
representations, as shown by perturbation theory up to 
two loops p2[ . Interest in Casimir scaling was recently 
revived it has been triggered by numerical studies 

of SU(3) lattice gauge theory [^|j4j, which indicate that 
Monte Carlo data for the potential between charges in 
different representions are consistent with Casimir scal- 
ing up to a relatively large distance, r « Ifm. 

The Casimir scaling law holds exactly in two- 
dimensional QCD. In higher dimensions no strong argu- 
ments exist in favor of a mechanism preserving Casimir 
scaling across the roughening transition, from strong 
to weak coupling; nor from small distance (essentially 
perturbative, characterized by a Coulombic potential) 
to large distance (characterized by a string tension for 
sources carrying charge). We have shown explic- 
itly that Casimir scaling does not survive the next- 
to-leading order calculation of the ratios R{k,N) in the 
strong-coupling lattice Hamiltonian approach. 

It is interesting to note that the sine formula (||) 
emerges also in the context of the two-dimensional 
SU(7V) X SU(7V) chiral models. As amply discussed in the 
literature (see e.g. Refs. ]l^ , p^ and references therein), 
d-dimensional chiral models and 2(i-dimensional lattice 
gauge theories manifest deep analogies in the continuum 
and on the lattice. Indeed, one may establish the follow- 
ing correspondence table for the lattice formulations: 



Chiral Models 

site, link 
loop 
length 
mass M 
two-point correlation 



Gauge Models 

link, plaquette 
surface 
area 
string tension a 
Wilson loop 



One may also add to this table the bound state masses 
Mk of chiral models and the fc-string tensions ak of gauge 
theories. In particular, in the case d = 1 the relation 
is exact, and one can prove that Casimir scaling holds 
for the masses of the bound states. In analogy to four- 
dimensional SU(iV) gauge theories, in two-dimensional 
SU(7V)x SU(A^) chiral models the Casimir scaling law 
holds for the strong-coupling limit of the corresponding 
lattice Hamiltonians (but it is not satisfied by the correc- 
tions) and for the small-distance behavior of the corre- 
lation functions related to different representations. On 
the other hand, the exact S-matrix [ p^ , derived using the 
existence of an infinite number of conservation laws and 
Bethe Ansatz methods, dictates that all bound states be- 
long to the rank-fc antisymmetric representations and sat- 
isfy the sine formula, Mk — Af sin(fc-^)/ sin(-^), where 
Alk is the mass of the fc-particle bound state. The 
question arises again: does this result extend to four- 
dimensional SU(-/V) gauge theories? 



This issue can be investigated numerically using the 
lattice formulation of SU(A^) gauge theories. Recent nu- 
merical results for R{2, N), obtained for = 4, 5 
show that i?(2, A^) < 2 ; thus, (T2 < 2cr, indicating that 
flux tubes attract each other. However, the error esti- 
mates on i?(2, A^) do not allow to exclude any of the 
two above-mentioned hypotheses. Indeed the sine and 
Casimir formulas give numerically close predictions for 
fc = 2, so that high accuracy is necessary to distinguish 
them. 

In this work we further investigate the spectrum of the 
string states. We present results from Monte Carlo (MC) 
simulations of the four-dimensional SU(6) lattice gauge 
theory using the Wilson formulation. For N = 6 there are 
two nontrivial fc-string tensions besides the fundamental 
one, thus providing a stringent test of the various ideas 
discussed above. We anticipate here our final results for 
the two independent /c-string tension ratios: 



i?(2,6) 1.72 ±0.03, 
i?(3,6) = 1.99 ±0.07. 



(4) 
(5) 



They are both consistent with the predictions of the sine 
formula (|), which are 5'(2,6) = 1.732... and S'(3,6) = 2, 
respectively. On the other hand, our results show devi- 
ations from Casimir scaling; the predictions in this case, 
C(2,6) = 1.6 and C(3,6) = 1.8, are off by approximately 
four and three error bars, respectively. 

TABLE I. MC results for the fc-string tensions. 
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lattice 




(j-ija 


0-3/0- 


0.342 


8^ X 20 


0.3151(6) 


1.65(2) 


1.91(3) 




12^ X 24 


0.3239(8) 


1.66(3) 


1.91(9) 


0.344 


12^ X 24 


0.2973(5) 


1.73(2) 


1.95(5) 


0.348 


10^ X 20 


0.2534(6) 


1.73(3) 


2.08(10) 




12^ X 24 


0.2535(6) 


1.71(4) 


2.06(11) 


0.350 


12^ X 24 


0.2380(6) 


1.72(3) 


1.95(9) 


0.354 


12^ X 24 


0.2103(5) 


1.73(3) 


2.04(6) 



In our simulations we employed the Cabibbo-Marinari 
algorithm ||l9| to upgrade SU(A^) matrices by updating 
their SU(2) subgroups (we selected 15 subgroups). This 
was done by alternating microcanonical over-relaxation 
and heat-bath steps, typically in a 4:1 ratio. Table | con- 
tains some information on our MC runs: The coupling 
values 7 = /3/(2A^^), lattice sizes, and the results 
for the /c-string tensions. The number of sweeps per run 
was typically above 500k, and measurements were taken 
every 10-20 sweeps. The values of 7 were chosen to lie be- 
yond the first order phase transition which occurs in the 
Wilson formulation of SU(A^) gauge theories for A^ suf- 
ficiently large (see e.g. Rcfs. [^,^). In order to deter- 
mine the value 7c where the first order transition occurs, 
we performed simulations starting from hot and cold con- 
figurations to display hysteresis, and from mixed-phase 
configurations, obtaining the estimate 7c = 0.3389(4). 
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We used asymmetric lattices (L^ x T) with a larger time 
size. For some values of 7 we performed simulations for 
two lattice sizes, to check for finite size effects. The lat- 
tice sizes L were chosen so that L^Ja > 2.5, and for most 
of them Ly/a « 3. This requirement ensures that finite 
size effects on fc-string ratios are negligible, as can be seen 
by comparison of the results for different sizes (see also 
Refs. jTzl). Further confirmation comes from preliminary 
results (« 200k sweeps) on a 16^ x 32 for 7 = 0.350. 

In our simulations we have also measured the topolog- 
ical charge Q, by a cooling technique. A severe form of 
critical slowing down is observed in this case: The auto- 
correlation time Tq for Q appears to increase exponen- 
tially, Tq cx exp(c/CT^/^) with c ~ 2.4. As a consequence, 
the run for the largest value of 7 considered, 7 = 0.354, 
did not correctly sample Q, presumably because it was 
not sufficiently long 300k sweeps). This dramatic ef- 
fect was not observed in the correlations used to deter- 
mine the fc-string tensions (a blocking analysis did not 
show significant time correlations in measurements taken 
every 10-20 sweeps), suggesting an approximate decou- 
pling between the topological and nontopological modes. 
This suggestion is also supported by the fact that string 
tension results for 7 = 0.354, extracted from a simulation 
which did not sample correctly Q, turn out to be in agree- 
ment with those for smaller 7, for which Q was sampled 
correctly. We note incidentally that this phenomenon has 
already been observed in simulations of two-dimensional 
CP^-i models |2|. 

The /c-string tensions are extracted from the large time 
behavior of correlators of strings in the antisymmetric 
representations, closed through periodic boundary con- 
ditions (see e.g. Refs. [ p3|Jl7| ): 

Cr{t) = (Xr[^(0; 0)] Xr[P{xuX2-M (6) 

where P(xi, a;2; i) = Iix3Uj,{xi,X2,x^]t), U{:x.;t) are the 
usual link variables, and Xr is the character of the rep- 
resentation r. In particular, x/[^'] = TrP for the funda- 
mental representation, and 

Xfe=2[P]=Trp2_(TrP)2, (7) 
Xk=3 [P] ^ 2Tr - 3Tr Tr P + (Tr Pf (8) 

for the antisymmetric representations of rank k — 2 and 
/c = 3, respectively. 

These correlators decay exponentially as exp(— mj-t) 
where is the mass of the lightest state in the cor- 
responding representation. For a /c-loop of size L, the 
/c-string tension is obtained using the relation ||23|| 

TT 

TOfe=0-fcP-— . (9) 

The last term in Eq. is conjectured to be a universal 
correction, and it is related to the universal critical be- 
havior of the flux excitations described by a free bosonic 
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FIG. 1. The scaling ratio i?(2, 6) as a function of a^a. 

string 1^,^. In order to improve the efficiency of the 
measurements we used smearing and blocking procedures 
(see e.g. Refs. |25)) to construct new operators with a 
better overlap with the lightest state. The masses irik 
were obtained by fitting the time dependence of the cor- 
relations. The fitting range is source of systematic error; 
we have checked that all reasonable choices of this range 
yield consistent results within the quoted errors. More 
details on the Monte Carlo simulations and the analysis 
of the data will be reported elsewhere p^ . 

The results for the ratios P(2, 6) and P(3, 6) are pre- 
sented in Table |, and plotted in Figs. |^ and || versus a^a 
to evidentiate possible scaling corrections, for the case 
k ^ 2 and k ~ 3 respectively, together with the sine for- 
mula (H) and the Casimir scaling predictions. The ratio 
R{2, 6) shows good scaling for 7 > 0.344. Scafing de- 
viations are observed only for 7 — 0.342, and this may 
be due to the vicinity of the phase transition. The data 
for 7 > 0.344 are consistent with a constant, thus we do 
not attempt to fit the dependence of our result on the 
lattice spacing a. Our final value for R{2, 6) is obtained 
by combining the results at 7 = 0.348 (for the largest 
lattice) and 7 = 0.350. The error we report is given by 
the typical error of each single point. Of course, this esti- 
mate assumes that the scaling corrections are small and 
negligible for a^a ~ 0.05. The data for smaller 7, and in 
particular the one for 7 = 0.344, are essentially used to 
check this fact. They suggest that the scaling corrections 
are at most of the same size of the error we report. MC 
runs at the largest value of 7, i.e. 7 = 0.354, show the 
aforementioned decoupling of the string tensions from the 
topological degrees of freedom; however, given the poor 
sampling of the topological charge in those runs, we do 
not include them in the final estimate of the string ten- 
sion ratio. Similar comments apply to the P(3, 6) ratio. 

We have also explored correlators in the symmet- 
ric rank-2 representation, finding no evidence for stable 
bound states, in accordance with general arguments and 
with the spectrum of chiral models. 

Our final estimates have been reported in Eqs. (0) and 
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FIG. 2. The scaling ratio R{3, 6) as a function of a^a. 



(g). They show deviations from Casimir scahng [g6[. It 
is worthwhile to emphasize that such corrections are to 
be expected, as mentioned previously. This fact is fur- 
ther confirmed by the computation of the ratio Ufc/cr to 
0{g~^) in the strong-coupling expansion of the lattice 
Hamiltonian formulation of d-dimensional SU(7V) gauge 
theories. We obtained lllsi 



Ok ^ k{N^k) 
a 



N -1 



1 



{d-2)f{k,N) 



(10) 



where /(fc, N) is explicitly fc-dependent. In particular 
/(2,iV) = f +0(^).^ 

In conclusion, we claim that our numerical results for 
the four-dimensional SU(6) gauge theory are consistent 
with the sine formula, and the universality hypothesis 
that is behind it. Of course, they do not prove that it 
holds exactly. But they put a stringent bound on the size 
of the possible corrections. On the other hand, our re- 
sults show a clear evidence of deviations from the Casimir 
scaling. This fact should be relevant for the recent de- 
bate on confiniment models, such as those discussed in 
Refs. However, Casimir scaling may still be con- 

sidered as a reasonable approximation, since the largest 
deviations we observed were about 10%. 

One last remark regards the large- TV behavior of the 
sine formula: S{k, oo) = k + O (l/A^^). In this respect 
the sine formula is peculiar because there are no a priori 
reasons for the /c-string tension ratio to be even in 1/N. 
The same observation applies to the two-dimensional 
SU(A^)x SU(7V) chiral models, but there we know that 
the sine formula holds and it comes from the structure 
of the S-matrix, which is essentially determined by the 
existence of an infinite number of conservation laws. 
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In order to further check the robustness of our statement 
concerning the deviation from the Casimir formula, we 
have analyzed our data using also fits which allow for 
scaling corrections. Of course, for this purpose one must 
consider all data, including the one at 7 = 0.342, since 
it is the only one showing apparent scaling violations. As 
suggested in Ref . jlTj , one may fit the data to the linear 
behavior A + Ba. The results are R{2, 6) = 1.77(5) and 
R{2, 6) — 1.79(6) respectively with and without the point 
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at the largest value of 7, i.e. 7 = 0.354. Therefore, if one 
wants to be more cautious in treating the systematic er- 
ror due to scaling corrections, one may take into account 
the hnear fits, and arrive at R{2, 6) = 1.72 ± 0.03 ; 
which covers all the results obtained above. This conser- 
vative estimate is still not consistent with the Casimir 
formula. Thus, our conclusions are fully justified even by 
this overly cautious analysis. 



